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A coherent control scheme for the population distribution in the vibrational states of nonpolar 
molecules is proposed. Our theoretical analysis and results of numerical simulations for the interac- 
tion of the hydrogen molecular ion in its electronic ground state with an infrared laser pulse reveal 
a selective two-photon transition between the vibrational states via a coupling with the first excited 
dissociative state. We demonstrate that for a given temporal intensity profile the population trans- 
fer between vibrational states, or a superposition of vibrational states, can be made complete for a 
single chirped pulse or a train of chirped pulses, which accounts for the accumulated phase difference 
due to the AC Stark effect. Effects of a spatial intensity (or, focal) averaging are discussed. 

PACS numbers: 33.80.Wz, 37.10. Mn, 42.50.Hz 



I. INTRODUCTION 

Control of a molecular process and, in particular, con- 
trol of the quantum state of a molecule is a long-standing 
goal in molecular physics, chemical physics, chemistry, 
and, more recently, strong-field physics and quantum in- 
formation technology. Light fields generated from fem- 
tosecond laser pulses have been tailored to manipulate 
the population in a molecule to a specific target state. 
Applications range from the control of the motion and 
shape of a wavefunction [1] to the control over molecu- 
lar fragments in a simple reaction '5]. Moreover, novel 
quantum technologies such as ultracold chemical reac- 
tions [3], frequency metrology of fundamental constants 
[3] or molecular quantum bits for quantum information 
processing [5] require the control over the internal quan- 
tum states of a molecule. In many of these applications 
it is desirable to completely evolve the molecular system 
in a robust way to the desired target quantum state. 

Much progress in the preparation and manipulation of 
the quantum state of neutral diatomic polar molecules [B]- 
[TT] and molecular ions OES] has been achieved recently. 
Several experimental approaches, e.g. photoassociation 
of ultracold atoms involving Feshbach resonances, spon- 
taneous decay and optical schemes such as stimulated 
Raman adiabatic passage (STIRAP), have been devel- 
oped. Here, we present a theoretical study and results of 
numerical simulations exploring an alternative method to 
manipulate and control the vibrational state population 
in a molecule. We performed our analysis for the sim- 
plest molecular ion, namely the hydrogen molecular ion. 
However, our results likely hold for other molecules as 
well. Since Hj is a nonpolar molecular ion, optical dipole 
transitions between vibrational states are forbidden. Fur- 
thermore, the first excited electronic state (2p(T„) is a 
dissociative state and higher electronic states are rather 
weakly bound. Thus, control and manipulation of the 
internal quantum state of H^ pose a challenge, since the 
application of methods such as STIRAP is limited due to 
the reduction in efficiency for a transition via continuum 




FIG. 1: (Color online) Coherent control scheme for pop- 
ulation transfer between vibrational states of a nonpolar 
molecule such as HJ. Vibrational levels (here, v — Q and 
1/ = 6) of the ground electronic state (Iso-g) of are cou- 
pled via multiphoton transitions with the first excited dis- 
sociative state (2p(T„). Efficient control over the population 
transfer between the states can be achieved for near-resonant 
net two-photon transitions. 



states, e.g. molecular dissociative states [H] . 

A sketch of the control scheme is shown in Fig. [T] We 
consider to transfer the population between two vibra- 
tional states {v = and v = as an example in Fig. [T]) 
in the electronic ground state of H^ [Isag). Both states 
are coupled by an intense infrared laser pulse via multi- 
photon transitions with the first excited dissociative state 
{2pau)- As we will show below, population is transfered 
from one vibrational state (or, a superposition of states) 
to the other if the photon energy approximately equals 
half of the energy difference between the two vibrational 
states. The transfer can be made complete by applica- 
tion of a chirped laser pulse or via coherent accumulation 
using a train of chirped pulses, which take account of the 
presence of the dynamic Stark shifts. 

The present scheme resembles analogy with well- 
known processes such as STIRAP |T5] , laser-induced con- 
tinuum structure (LICS) [16j or Raman chirped adiabatic 
passage (RCAP) |T7j but it also differs in a number of 
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characteristic features. For example, the energy differ- 
ence between the excited state and the initial or final 
states is much larger than the photon energy and, con- 
sequently, the transition via the excited state is a mul- 
tiphoton transition. Traditionally, these kind of coher- 
ent control schemes are based on one-photon transitions 
between the states involved. Related to this difference, 
application of the present scheme is based on the interac- 
tion with just one laser pulse, while often two laser pulses 
operating at different frequencies are used. As in the pre- 
vious schemes, maximum population transfer is achieved 
for a certain laser frequency. Interestingly however, the 
related condition corresponds to a two-photon resonance 
between the initial and the final state. 

The rest of the paper is organized as follows: We first 
present in Section II. A. the theoretical model for the hy- 
drogen molecular ion, which we used for the numerical 
simulations of the population transfer between two vi- 
brational states. Based on the numerical results we will 
identify the condition needed to maximize the population 
transfer in the field of an unchirped pulse, namely a slight 
detuning from the two-photon resonance between the ini- 
tial and final states. Next, in section II. B., we analyze 
the process using a simplified model which takes account 
of just two electronic states. The analysis shows that the 
transfer can be made complete using a chirped pulse and 
an analytical formula for the chirp will be derived. In 
section II. C. we verify our analytical predictions by per- 
forming calculations using the two-state model and the 
full numerical simulations for transfers between individ- 
ual states as well as superposition of states. Results of 
both set of calculations are in perfect agreement show- 
ing the reliability for our analysis based on the two-state 
model. Finally, we investigate the effect of a spatial in- 
tensity (or, focal) averaging before we conclude with the 
main findings presented in this paper. 



II. THEORETICAL MODELS AND RESULTS 

We analyzed the control scheme using analytical 
calculations as well as numerical simulations of the 
time-dependent Schrodinger equation (TDSE). Wc first 
present the results of our numerical simulations for inter- 
action of the molecular ion with an unchirped pulse at 
wavelengths in the infrared. This will lead us to an anal- 
ysis of the population transfer within a two-state model. 
Based on the conclusions of our model analysis we will 
then verify by performing calculations using both the 
two-state model as well as the full ab-initio numerical 
simulations that the transfer can be made complete for 
a chirped pulse. 



A. Two-photon resonance condition 

For the numerical simulations we used a two- 
dimensional model of the hydrogen molecular ion, in 



Vibrational State 


Energy (a.u.) 


i/ = 


-0.776 


V = 1 


-0.767 




-0.758 




-0.749 




-0.741 



TABLE L Energies of the first vibrational levels in the ground 
electronic state of the model used in the present study. 



which the motion of the electron as well as the motion 
of the protons are restricted along the polarization di- 
rection of the external field. The Hamiltonian of the 
system is given in length gauge as (Hartree atomic units, 
e — m = h = 1, are used throughout) [T5] : 



H{x,R,t) 
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where (P, R), and (p, z) are momentum operators and 
positions of the relative coordinate of the two protons 
and the electron, respectively. /Zp = A//2 is the reduced 
mass with M — 1836 is the mass of the proton, and 
ckg = 1 and Up ~ 0.03 are soft-core parameters |18| . 
E{t) = Eq cos{ujt) SIT? {irt / NT) is the electric field of the 
linearly polarized laser pulse. The corresponding TDSE 
was solved using the Crank-Nicolson algorithm on a grid 
with spatial spacings of Ai?=0.03 au and Az^O.l an, 
the time step was less than At= 0.05 au. The vibra- 
tional levels in the ground electronic state of the model 
were found by imaginary time propagation and the cor- 
responding energies are given in Table 1. 

First, we performed a series of calculations in which 
we prepared the model system in the ground vibra- 
tional state. We tried to find the optimal wavelength of 
the light field for a selective population transfer into one 
of the excited vibrational states. In Fig. [2] we present 
our results for the final populations in the first four ex- 
cited states driven by laser pulses at A = 9919.9 nm, A 
= 5059.3 nm and A = 3441.0 nm. The photon energy is 
equal to half of the energy difference between the first, 
second, and third excited state and the ground state, re- 
spectively. The pulse length was 10 field cycles and the 
peak intensity was (a) 10"^^ W/cm^, and (b) 10"^^ W/cm^. 

The results show that a particular vibrational state can 
be selectively excited via a (net) two-photon transition 
from the initially prepared ground state. The popula- 
tions in all the other states are at least three orders of 
magnitude smaller as compared to the maximum excited 
state population. However, we observe a transfer of pop- 
ulation of about 1% only, or even less. We note that for 
long pulses we observe Rabi-like oscillations but without 
surpassing the 10% level for the population in the excited 
state. 
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FIG. 2: (Color online) Final population of the first vibra- 
tional states by interaction with laser pulses at 9919.9 nm, 
5059.3 nm and 3441.0 nm, respectively. Results of the pro- 
jections to = 5 and u = 6 are below 10~^^ and not shown. 
The number of cycles was 10 in each simulation, and the peak 
intensity was (a) 10^^ W/cm^ and (b) 10^^ W/cm^ 



Next, we present in Fig. [3] results for the population in 
the V — \ state as a function of the laser wavelength for 
two different peak intensities, (a) 10^^ W/cm^ and (b) 
10^^ W/cm^. The number of cycles is kept the same as 
in Fig. [2j The results show that the population transfer 
is maximum for a slight detuning A from the two-photon 
resonance condition wi^o = {E„=i — -Ey=o)/2. The de- 
tuning is found to depend linearly on the intensity (see 
Fig. |3jc)). These findings appear to hold in general in- 
dependent of the final state, since we find the same linear 
dependence for excitation of the 1^ — 2 state (Fig. |4| . 



B. Analytical model 

Based on the findings from our numerical simulations 
we would like to address the following questions: What 
is the physical mechanism of this population transfer be- 
tween vibrational states in a nonpolar molecule? What 
are the origins of the near-resonant two-photon condi- 
tion and the detuning? And, finally, is it possible to 
enhance the efficiency of the population transfer? In or- 
der to answer these questions, we developed the following 
analytical model. 

We assume that the photon energy is small as com- 
pared to the energy difference between the electronic 
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FIG. 3: (Color online) Final population of i/ — 1-state as 
a function of the laser frequency. Note that the maximum 
occurs for a slight detuning A from the two-photon resonance 
frequency ujifi = (Ev=\ — Ev^o)/2 au. The number of cycles 
was 10 in each simulation and the peak intensity was (a) 10^^ 
W/cm^ and (b) 10^^ W/cm^. (c) Detuning A as a function 
of intensity. 
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FIG. 4: (Color online) Same as Fig. |3j but for u = 2-state. 

states and the peak intensity is sufficiently small such 
that the state of the molecular ion, which is initially pre- 
pared in a superposition of vibrational states of the elec- 
tronic ground state (IsCTg), is well described by a super- 
position of the two lowest electronic states at any time 
t during the interaction (in Born-Oppenheimer approxi- 
mation): 

+0,„(x,i?) / dkbk{t)^Jk{R) ■ (2) 



This reduction to a two-level model is justified within the 
range of laser parameters considered here, since we do 
not observe any ionization, dissociation or population to 
higher excited states in our full numerical simulations. As 
we will show in section II. C below, the results obtained 
using the two-level model do agree very well with those 
of the full numerical simulations. 

Inserting Eq. ^ in the Schrodinger equation 
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id'^{x,R,t)/dt = H{x,R,t)^{x,R,t), multiplying by 
<l)%^{x,K)i)l{R) or (l)%^{x,R)^l{R) from the left and in- 
tegrating with respect to x and R we obtain a set of 
equations for the amplitudes au{t) and bk[t): 



ibk{t) = E,^,kbk(t)+£{t)Y, f^t,kaAt)' 



(4) 



with ^i^^k = {4',yg{x,R)i)^{R)\x\(j)c,S^,R)ipk{R)), where 
the brackets indicate integrations over x and R. £'0-3.1/ 
and £'cr„,fc are the field- free energies of the states. Eqs. 
([3]) and Q represent the coupling of the vibrational 
states in the electronic ground state with the dissocia- 
tive state. Since we are interested in the evolution of the 
population in the vibrational states, we integrate Eq. Q 
with respect to time and then substitute the result into 
Eq. ^ to get: 



dk nl, n^^k X 



dt'£{t)£{t')c^,{t')e 



iAE„ 



a„(t)e'-^"9'-* with |a„(t)| 



(5) 

Cv{t)\ and 



where c^{t) 

AEy^k — £cr„,fc ~ ■^<^g,v- If the laser frequency is far 
off-resonant, i.e. /S.E^,k 3> uj, we can assume that the 
electric field strength £{t) and c^(i) are slowly varying in 
time compared to the exponential term, which depends 
on the energy difference l^E^^k, and perform the integral 
in Eq. ([5]) approximately to get: 

where /ij^, = J dk fi^, f. fi^^k/AE,y'^k- For the sake of 
simplicity we defined E^^^i, = E^. The temporal evo- 
lution of the amplitude Ci, depends quadratically on the 
electric field strengths, which gives rise to the two-photon 
transition between the vibrational states observed in the 
numerical simulations. 

Please note that the current control mechanism is 
based on muUiphoton transitions for which the off- 
resonant condition AEy^k 3> io holds. This is the main 
difference between the present mechanism to induce tran- 
sitions between different internal states of the molecule 
and processes such as LICS (laser-induced continuum 
structure, [16]) and STIRAP [15]. For example, in STI- 
RAP the off-resonant condition AEy j^ ^ w does not hold 
and Eq. ^ cannot be further simplified as in the present 
analysis. 

Using the two-level approximation between an initial 
(i^i) and a final (i/y) vibrational state we obtain: 




JAEt 



2 -lAEt 




•(7) 



where AE = E^^ - E^^. Note that [i^^^^j ^ Ait-^.i^.i if 
AE^.^k 7^ AE^j^k. In general, this gives rise to a non- 
Hcrmitian Hamiltonian in Eq. ([T]) and, hence, a loss 
of population in the two-level system via excitation into 
the dissociative state. The loss becomes negligible, if 



ia^{t) = E^^^^a^{t)+E{t) dk ^JL^^kbk(t), (3) AE ^ AE, 



i,._k ~ AE^j._k- This is the case for the (low) 
vibrational states in H J , correspondingly we did not ob- 
serve any dissociation of in our numerical simulations. 

Transforming c„^{t) = exp{i J^^£^{t')fil^,^Jt')c^^{t) 
(with j = i, /) yields: 













(8) 



with 



(9) 

(with j — i, /) represents the time-dependent Stark 
shift of the vibrational state j, which results in a phase 
shift between the states during the interaction with the 
field. A similar result was obtained in Ref. [1^1 HO] 
for transitions between atomic bound electronic states. 
Please also note that for nl.^^. — Eq. (|9| re- 

duces to the standard two-level system coupled via a 
two-photon transitions with an electric field [21j . 

According to the analysis above we expect that a com- 
plete coherent transfer of the population from one vibra- 
tional state to the other can be achieved by designing a 
chirped pulse in which the frequency variation follows the 
time-dependent accumulation of the Stark phase shift as: 



2 



2u{t)t = AEt 

= 2{ujf,,t + 5{t)) , 



2 



£^{t')dt' 



(10) 



where ujf^i satisfies the two-photon resonance condition. 

In general „ — v' unknown, in our numerical 
model as well as in an experiment. However, by using 
a chirped pulse with 5{t) = a J*^£^{t')dt' and sweeping 
the constant a one can optimize the transition probabil- 
ity. We may note parenthetically that this method may 
be also applied to gain information about an unknown 
molecular system, which we however will not further ex- 
plore here. 



C. Population transfer with chirped pulses 

In this subsection we will test our prediction that the 
population transfer can be made complete for a certain 
chirp of the pulse, using both the two-level model as 
well as full numerical simulations. To this end, we de- 
termined the chirp parameters for the transition from 
the vibrational ground state (f^ = 0) to the second ex- 
cited state {vf = 2) as — ^2.2 = (—2.66 ± 0.02) 
au and 2 = (0.255 ± 0.001) au, and performed full 
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FIG. 5: (Color online) Population in the v = 2 state during 
the interaction with a chirped laser pulse (280 cycles, dashed 
line) and a train of 28 chirped pulses (10 cycles each, separated 
by 1 cycle, solid line), by solving (a) two-level matrix equation 
Q and (b) the TDSE. The peak intensity of each pulse was 
10" W/cm^ 



9919 nm 




FIG. 6: (Color online) Spatial intensity profile of a Gaussian 
laser pulse (scaled in 10^^ W/cm^). TDSE simulations for 
population transfer to (a-c) i/ = 1 using a chirped laser pulse 
with frequency oji^o (9919 nm) and 120 cycles and (d-f) u = 2 
using a chirped laser pulse with frequency uj^fi (5059 nm) 
and 280 cycles were performed. The results presented in the 
different panels correspond to the peak intensities labeled on 
the intensity profile: (a, d) lO" W/cm^ (b, e) 6.0 x 10^^ 
W/cm^ and (c, f) 1.3 x 10^^ W/cm^. The chirp parameters 
were determined to maximize the population transfer for the 
peak intensity at the center of the pulse and kept the same 
for all simulations, wo refers to the beam waist. 



time-dependent numerical simulations as well as calcula- 
tions using the two-level approximation for our model 
system. The results of both approaches agree very well 
(dashed lines in Fig. [5]). The final probability to find the 
system in the excited vibrational state is indeed almost 
100% (exactly, 99.89%). Since the excitation probabil- 
ity depends on the square area of the pulse envelope, 
the same result can be achieved by a coherent accumula- 
tion of the transition probability using a train of chirped 
pulses (solid lines in Fig. [5|. The delay between subse- 
quent pulses must be such that the phase of the initial 
state at the beginning of a pulse equals the phase at the 
end of the previous pulse [55] . 

Please note that our condition for the chirp parameters 
depends on the temporal intensity profile, in particular 



FIG. 7: (Color online) Population transfer from a superpo- 
sition of the excited vibrational states !/ = 1 and u = 2 to 
the vibrational ground state, (a) Population of the lowest 
vibrational states during the interaction with a first chirped 
pulse (the same used in Fig. [H] ijJ2,o = 9.00 x 10~^). (b) A 
second chirped pulse with a peak intensity of 10^^ W/cm^, 60 
cycles and /xq.q — — (—1.56 ± 0.02) au interacts with the 
remaining superposition after a time delay 10.8 fs. 



the maximum intensity (c.f. Eq. (lOl). In an exper- 
iment averaging over the spatial intensity profile (also 
called focal averaging) cannot be avoided. This results 
in a variation of the maximum intensity in the temporal 
profiles over the focal area. In order to analyze the effect, 
we have considered a Gaussian (spatial) intensity profile 
with a maximum peak intensity of 10^^ W/cm^ (see. Fig. 
|6]). We determined the chirp parameters for a complete 
population transfer from the vibrational ground state to 
the first (left hand panels in Fig. [6| and second vibra- 
tional state (right hand panels in Fig. |6| for the peak 
intensity at the center and performed TDSE simulations 
for a set of peak intensities keeping the chirp parameters 
the same. As expected, the population transfer at the 
center is complete but it becomes less and less efficient 
as the intensity drops towards the wings of the intensity 
distribution. Still, we find a transfer of more than 73% 
(for the transition to the first excited state) and more 
than 50% of the population (for transfer in v = 2) at 
60% of the maximum peak intensity. Even, in the wings 
of the intensity distribution the transfer is with more of 
1%, not negligible. A similar reduction in the control effi- 
ciency will occur for an unaligned ensemble of molecules, 
since the effective interaction strength does scale with the 
cosine of the angle between the polarization direction of 
the field and the internuclear axis. 



Finally, we considered population transfer from a su- 
perposition of the first and second excited vibrational 
level {i> = 1 and i/ = 2) to the vibrational ground 
state. The results of our numerical simulations, shown 
in Fig. [T] demonstrate that with two chirped pulses with 
W2,o = 9.00 x 10-3 au (5059.3 nm) and wi^o = 4.59 x lO^^ 
au (9919.9 nm) and the appropriate chirps first the pop- 
ulation from the second vibrational state (Fig. [Tf^a)) and 
then the remaining population from the first vibrational 
state (Fig. [Tjjb)) is transferred to the ground state. 
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III. CONCLUSIONS 

In conclusion, we have theoretically and numerically 
analyzed the coherent control of the population distri- 
bution in the vibrational states of the nonpolar hydro- 
gen molecular ion. Our results unveil that a selective 
excitation of vibrational states can be achieved via a 
(net) two-photon transition using infrared laser pulses. 
The transitions are accompanied by a dynamical Stark 
shift, which depends linearly on the laser intensity. Us- 
ing a chirped laser pulse the Stark phase shift can be 
followed and for a given temporal intensity profile a com- 
plete coherent transfer can be achieved between two vi- 
brational states of the molecule. Looking ahead, we an- 
ticipate that the same coherent control scheme can be 
applied to other nonpolar molecules with bound or dis- 
sociative electronic states as well as to the preparation of 
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